Low-energy models for correlated materials: bandwidth renormalization from Coulombic screening 
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We provide a prescription for constructing Hamiltonians representing the low energy physics of correlated 
electron materials with dynamically screened Coulomb interactions. The key feature is a renormalization of the 
hopping and hybridization parameters by the processes that lead to the dynamical screening. The renormaliza- 
tion is shown to be non-negligible for various classes of correlated electron materials. The bandwidth reduction 
effect is necessary for connecting models to materials behavior and for making quantitative predictions for 
low-energy properties of solids. 
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A key step in the theoretical analysis of strongly correlated 
materials is the derivation, from an all-electron Hamiltonian in 
the continuum, of an effective model which correctly captures 
the physics of the low-energy degrees of freedom. Tremen- 
dous progress in this direction has been achieved by using 
density functional theory (DFT) techniques [Q to compute a 
full set of energy bands, from which a subset of correlated 
orbitals is abstracted for further detailed study using many- 
body (typically dynamical mean field (DMFT)) methods ||2}- 
4 1 or LDA+U. The interaction parameters used in the many- 
body studies are the matrix elements of the screened Coulomb 
interaction in the correlated subspace. Various methods are 
used to obtain the screened matrix elements, including the 
constrained local density approximation [5], linear response 
[|6|, or the constrained random phase approximation (cRPA) 
Q. This DFT+cRPA+DMFT approach enables quantitative, 
testable theoretical predictions for correlated materials. 

In this paper we show that this scheme misses an impor- 
tant aspect of the physics: the downfolding produces a dy- 
namically screened Coulomb interaction which leads to an 
effective model with a bandwidth that is reduced relative to 
the starting (e.g. DFT) bandwidth and a low energy spectral 
weight which is also reduced. This effect has previously been 
noticed QE). A similar renormalization was also discussed 
in the context of Holstein-Hubbard models in Refs ll9l [lOl . 
We present an explicit nonperturbative prescription for deter- 
mining the renormalizations quantitatively, and demonstrate 
that the resulting effective model provides a good descrip- 
tion of the low-energy part of the full (dynamically interact- 
ing) model over wide parameter ranges. Computations of the 
renormalizations for wide classes of correlated electron mate- 
rials indicate that their inclusion is crucial for a quantitative 
description, in particular resolving a long-standing discrep- 
ancy between the cRPA estimate of the Coulomb interaction 
and the values needed to describe experiments. 

We first provide a demonstration for the simplest case, 
where the downfolding from the full band structure is to 
a one-band model with hopping amplitude tij between the 
lattice sites i and j. Electrons with spin a in the cor- 
related orbital localized at site i are created [annihilated] 



by the operator d ia . [di a ]. Double occupation of a given 
atomic site costs a Coulomb energy U, which is renormal- 
ized from a bare value V (obtained from the site-local ma- 
trix elements of e 2 jr among the correlated orbitals) because 
of screening by degrees of freedom eliminated in the down- 
folding process. The interaction thus takes the general form 
\ (V5(t) + f/ ret (T)) n i (r)n i (r'). Screening is contained in 
the retarded part U Iet , which is parametrized by a contin- 
uum of modes of energy v with coupling strength \ 2 {v) — 
— lmU Ki {v)/iT, determined by the charge fluctuations, 

/>oo 

U m (r) = - / d^A 2 (^)cosh[(r-/3/2)i/]/sinh[^/2]. (1) 
Jo 

For simplicity of presentation we assume at first that there is 
only one important bosonic mode of energy ujq and coupling 
strength A. The Hamiltonian is then 

ijcr i ia 

+ u Y / blb l + \Y / 4*diAb l + b\). (2) 

i ia 

A Lang-Firsov (LF) transformation ifTTI fT2l H ->■ Hr , F = 
e s He~ s with S = - X) i(7 n ia (fy + b\) allows one to 
rewrite the model in terms of the polaron operators c| ff = 
exp(A(6t _ bi ))dl and c ia = exp(^(6, - b\))d la . We 
note that c and c) obey the same fermionic anti-commutation 
relations as the original electronic operators (d and d^). Ne- 
glecting one-body terms which can be absorbed in a chemical 
potential shift, we have 

ij(T % 1 

with the screened Hubbard interaction Ur, = V — — . 

We now assert that the low energy effective model is given 
by the projection of Eq. Q onto the subspace of zero-boson 
states, H e ff — (0\H\Q). The effective model is then 

#eff = Z B Ujdl a d j(T + Uq 22 d if d it d il d il, ( 4 ) 
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that is, an effective Hubbard model with an instantaneous 
interaction corresponding to the low frequency limit of the 
screened interaction and a new feature, namely a bandwidth 
renormalized by Zb = cxp(— A 2 /wq). An additional physi- 
cal consequence of the low-energy projection is that the pho- 
toemission spectral weight in the frequency range described 
by the effective model is reduced by the factor Zb relative 
to what would naively follow from Hlf- Mathematically, 
glow-energy^ tne physical electron Green function in the fre- 
quency range described by the effective model, is 



half filling with Z B =0.741 



G 



low-energy 



(r) = -Z B (Td i (T)d](0)} Ba 



(5) 



where -(Tdi(T)d](0)) ffcfr is the Green's function G-*(t) of 
the effective Hamiltonian H e ff in Eq. (j4j). Thus the observ- 
able spectral function ^4 low - ener gy = -IlmG low ' energy (a; - iS) 
becomes 



^low-energy / ^ 



— ImG eff (a;-^). 

7T 



(6) 



The physical origin is that part of the physical photoemission 
spectrum corresponds to the simultaneous creation of a hole 
and a plasmon excitation; these plasmon shakeoff processes 
account for the remaining 1 — Zb spectral weight. 

The effective model becomes an exact description of the 
low energy physics only when the ratio of the boson frequency 
ujq to a relevant energy E* diverges, but we find that the ef- 
fective model gives a reasonably good description even for 
loq/E* not too large. As an example, Table [i] compares exact 
results (obtained using the methods of Ref. B131 ) for the criti- 
cal interaction strength C/ cr ; t needed to drive a metal-insulator 
transition in single-site DMFT to the predictions of the ef- 
fective model. In these computations we assume that the 
Uj give a semicircular density of states with half-bandwidth 
D = 1. Combining previously computed single site DMFT 
results [13] with our bandwidth reduction prescription gives, 
at inverse temperature (5 — 100/D, an effective model predic- 
tion Uf it fa 2.55Z B . One sees that the effective model result 
is within 15% of the exact result except when there is strong 
screening and the boson frequencies are smaller than the full 
bandwidth (2 in present units). 

Figurefllcompares the electron spectral function, calculated 
from Eq.fpl) with semicircular density of states (half band- 
width D = 1), for screened interaction Uq — 2 with values of 



TABLE I: Critical interaction strength U^" (presented in terms of 
zero frequency screened value) needed to drive the metal insula- 
tor transition obtained from the single-site DMFT approximation to 
Eq. |2| at inverse temperature /3 = 100 and compared to the es- 
timate U^ t for different values of the screening frequency ujq and 
strength A. Also shown is the Lang-Firsov renormalization factor 
Z B = exp[-A 2 /^]. 
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FIG. 1: Spectral functions computed from Eq. (j2} at various screen- 
ing frequencies ujq with j3 = 40, screened interaction Uo = 2 and 
coupling constants chosen to produce the renormalization factor Zb 
as indicated. Also shown are the spectral functions computed from 
the effective model (Eq. l|6]l) and for the static U approximation. 



Zb representative of typical correlated electron materials to 
two approximations: the effective model defined above, and a 
"static U model" which uses the static value of the screened 
Coulomb interaction but does not include the bandwidth re- 
duction. The static U model corresponds to what is normally 
done in DFT+DMFT calculations. The analytic continuations 
are obtained using the technique proposed in Ref. fl41 . We 
see that the effective model with bandwidth reduction Zb re- 
produces very well the effective bandwidths of the Hubbard 
bands for all w taken into account here, which vary from 10 
down to 2.5. Even the smallest ujq, which is not in the antia- 
diabatic regime, yields Hubbard bands qualitatively well de- 
scribed by the static model with bandwidth renormalization 
Zb- The static U model is seen to be a poor approximation. 



TABLE II: Quasiparticle residue a = 1/(1 — dlm[E(icj)]/duj\ UJ= o) 
computed from the effective Hamiltonian Eq. Q with screened Uq = 
2, and different ujq, Zb and particle density as shown. The values in 
parenthesis give the relative discrepancy \a(uJo) / a(ujQ — > oo) — 1|. 
Note that the static model without bandwidth reduction (last row) is 
substantially incorrect. 
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Table [II] shows the results of an alternative analysis, car- 
ried out at the level of the quasiparticle renormalization a — 
1/(1 — dT,/diuj n ), which is obtained directly from the imagi- 
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nary time computations. We see that the "static U" result gives 
renormalization factors in error by factors of two or more in 
the half filled, strongly correlated case, and also unacceptably 
large errors in the weakly correlated quarter filled case. The 
effective model (row ujq = oo) is very close to the exact result 
for all screening frequencies in the weakly correlated quarter- 
filled case and is reasonably close to the exact result even as 
the adiabatic limit is approached. 

Analogous arguments for a model comprising also itinerant 
p states, and thus hopping parameters T vv , T v d and Tdd lead 
to a renormalization of each d operator by the factor V 'Zb = 
(0| exp(^(&i — 6^))|0) so that the hopping part of the one- 
particle Hamiltonian is renormalized as 



(ptrft) 




ZBTdd J 



(7) 



where the site dependence of each orbital species is not explic- 
itly denoted. Equation (|7]i shows that the bandwidth reduction 
implied by our effective model happens in a non trivial way in 
the case of the multi-band models usually dealt with in first- 
principles calculations. 

The arguments we have given are readily generalized to the 
case of an arbitrary dynamical interaction. The representation 
of Eq. ([T]l corresponds to a continuum of boson excitations, 
bi(u), one for each frequency v in the screening process, with 
coupling \(y). We then apply a generalized LF transforma- 
tion obtaining 



/>oo 

U = V + 2/n dulmU Kt (v)/v, 
Jo 

Zb = cxp J dv lmU Ki {v) / 



(8) 
(9) 



Matching this to the single mode formula implies a character- 
istic frequency 



w = 



f^dis vlmU Kt (v)/v 2 
f^dv \mU K i{v)/v 2 



(10) 



Our theory has important implications for electronic struc- 



ture calculations for correlated materials. Table (IIIi presents 
our results for ojq, Zb and U values for a range of compounds 
calculated using the cRPA method [7 |, in the implementation 
of Ref. 15 Typical Zb values for oxides or pnictides lie in the 
range of ~ 0.6 — 0.7 indicating substantial renormalization 
of the low energy bandwidths relative to DFT calculations, 
even though the screening frequencies ujq are typically high. 
Standard DFT+DMFT calculations are available for most of 



the compounds. As shown in Table III in these calculations, 
obtaining agreement with experimental results for mass en- 
hancements and metal-insulator phase diagrams has required 
the use of U values substantially (~ 40%) larger than the low- 
frequency Hubbard interactions calculated from cRPA. For 
example, for SrV03, LDA+DMFT calculations with U rang- 
ing from 4 to 5 eV were found to yield good agreement with 
experiments lfT744T9ll (instead of the cRPA value of 3.5 eV). 
Similarly, in VO2, U — 4.0 eV was used Il24ll25l instead of 
U = 2.7 eV. We believe that the difference arises because the 



TABLE III: Boson renormalisation factor Zb, characteristic fre- 
quency cjo [eV], bare interaction V [eV], zero-frequency screened 
interaction Uo [eV] as calculated within the cRPA, in the implemen- 
tation of Ref. 1151 . For the oxide and sulfide compounds (except 
SrMnOa), data refer to a model comprising only the t2 S states, where 
U is defined as the average over the diagonal entries of the Hubbard 
interaction matrix U mm ,mm- For the pnictide compounds, as well as 
for SrMnOs and CuO, a hybrid "d-dp" model in the notation of Ref. 
1 1511 161 was constructed and U(= Fq) is defined as the average over 
all density-density interaction matrix elements. Experimental lattice 
structures (rutile in the case of VO2, hexagonal lattice in the case of 
TaS2) were used except in the cases of Sr2 VO4, LaVC>3 and SrMnC>3, 
where an undistorted (double) perovskite structure was adopted. The 
column headed C/ut gives U values obtained via a variety of methods 
other than cRPA claimed in the literature to give quantitative agree- 
ment with experiment when used in DFT+DMFT (oxides, sulfides 
and pnictides) or DFT+U calculations (SrMnC>3 and CuO) within 
the same correlated subspace, but without the band renormalization 
physics. 
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previous literature did not incorporate the bandwidth reduc- 
tion effect, and artificially compensated this by increasing U. 
The one apparent exception is SrMnC>3, where the U value 
quoted in Ref. |23]was chosen to be consistent with the mag- 
netic moment but gaps or other dynamical properties were not 
studied. A more recent work of a t2 9 -only model required a 
rather larger value of 3.5 eV, but overlap with e g bands pre- 
cludes a cRPA estimate of Zb in this case. 

Figure [2] shows another illustration of the bandwidth renor- 
malization phenomenon, comparing the spectral function of 
optimally doped BaFe2As2 obtained with the "static U" ap- 
proximation (panel (a)) to the full treatment of the dynamic 
U, as explained in Ref. 3 1 (panel (b)), and the effective model 
(panel (c)). Comparison of panels (a) and (b) shows that 
screening has a substantial effect on the band structure, shift- 
ing the energy positions of bands and band crossings to a sig- 
nificant extent. (The model with screening also has an in- 
creased broadening resulting from a change in proximity to a 
spin freezing line whose position depends very sensitively on 
parameters 1311 : this effect is not of primary interest here). 
Comparison of panels (b) and (c) shows that the effective 
model captures the changes in band energies very well, and 
also reproduces the change in lifetime. 

To summarize, in this Letter we showed that the low en- 
ergy effective Hamiltonian relevant to correlated electron ma- 
terials involves two renormalizations: a reduction, to a value 
smaller than the isolated atom value, of the on-site Coulomb 
interaction and a reduction, to a value smaller than the band 
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FIG. 2: fc-resolved spectral function for K I Ba( 1 _ I )Fe2As2 at optimal doping x = 0.4 and /3 = 20eV _ , reported for a static U standard 
DMFT calculation (panel (a)), the DMFT calculation with dynamic U (to) (panel (b)), and the DMFT calculation for our effective low-energy 
model. In all calculations, the static limit of U(= Fa) is U(0) = 2.84 eV, and J = 0.68 eV. In the effective model, the double counting 
correction is set to match the d-electron number of the dynamical calculation. 



theory value, of the bandwidth. The reduction of the on- 
site Coulomb interaction is a straightforward consequence of 
screening by high energy degrees of freedom and has been 
discussed in many works. The reduction of the bandwidth is 
a more subtle effect, which has important consequences for 
the low energy physics, including a reduction in the ampli- 
tude, and a narrowing of the width of the low-energy part of 
the electron spectral function, as well as a shift in the loca- 
tion of the metal-insulator transition. We have provided a pre- 
cise prescription for obtaining the bandwidth reduction and 
have tested our low-energy effective description against nu- 
merically exact dynamical mean field solutions of Hubbard 
models with full dynamic U in a range of parameters rele- 
vant for correlated materials. Important open questions are 
the issues of full charge-self consistency and the related dou- 
ble counting correction, both of which require knowledge of 
physics at energy scales above the range of validity of the low- 
energy effective model. This is the subject of current research. 

We acknowledge useful discussions with T. Miyake, H. 
Jiang and G. Sawatzky. This work was supported by the 
French ANR under project SURMOTT, IDRIS/GENCI un- 
der Grants 2012096493 and 201201393, SNF Grant PP0022- 
118866, FOR 1346 and the US DOE under grant BES ER- 
046169. 



[1] W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

[2] S. Biermann "LDA+DMFT" - a Tool for Investigating the Elec- 
tronic Structure of Materials with Strong Electronic Coulomb 
Correlations, in Encyclopedia of Materials: Science and Tech- 
nology. Ed. by J. Buschow, Elsevier 2005. 

[3] K Held et al. J. Phys.: Cond. Matter 20, 064202 (2008). 

[4] G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Par- 
collet, and C. A. Marianetti, Rev. Mod. Phys. 78, 865 (2006). 

[5] O. Gunnarsson, O. K. Andersen, O. Jepsen, and J. Zaanen, 
Phys. Rev. B 39, 1708 (1989). 



[6] M. Cococcioni and S. de Gironcoli, Phys. Rev. B 71, 035105 
(2005). 

[7] F. Aryasetiawan, M. Imada, A. Georges, G. Kotliar, S. Bier- 
mann, and A. I. Lichtenstein, Phys. Rev. B70, 195104 (2004). 

[8] Y. Imai, I. Solovyev, and M. Imada, Phys. Rev. Lett. 95, 176405 
(2005). 

[9] Y. Takada and A. Chatterjee, Phys. Rev. B 67, 081102 (R) 

(2003) . 

[10] A. Macridin, G. A. Sawatzky, M. Jarrell Phys. Rev. B 69, 
245111 (2003). 

[11] I. G. Lang and Y. A. Firsov, Sov. Phys. JETP 16, 1301 (1962). 
[12] P. Werner and A. J. Millis, Phys. Rev. Lett. 99, 146404 (2007). 
[13] P. Werner and A. J. Millis, Phys. Rev. Lett. 104, 146401 (2010). 
[14] M. Casula, A. Rubtsov, and S. Biermann, Phys. Rev. B 85, 
035115 (2012). 

[15] L. Vaugier, PhD thesis, Ecole Polytechnique (201 1), L. Vaugier, 

H. Jiang, S. Biermann, in preparation. 
[16] M. Aichhorn, et al., Phys. Rev. B 80, 085101 (2009). 
[17] A. Liebsch, Phys. Rev. Lett. 90, 096401 (2003) 
[18] E. Pavarini, S. Biermann, A. Poteryaev, A. I. Lichtenstein, A. 

Georges, and O. K. Andersen, Phys. Rev. Lett. 92, 176403 

(2004) . 

[19] A. Sekiyama, et al., Phys. Rev. Lett. 93, 156402 (2004). 

[20] F. Lechermann, A. Georges, A. Poteryaev, S. Biermann, M. 

Posternak, A. Yamasaki, and O. K. Andersen, Phys. Rev. B 74, 

125120(2006). 

[21] R. Arita et al., J. Phys. : Cond. Matter 19, 365204 (2007). 
[22] M. De Raychaudhury, E. Pavarini, and O. K. Andersen, Phys. 

Rev. Lett. 99, 126402 (2007). 
[23] J. Hee Lee and K. M. Rabe, Phys. Rev. Lett. 104, 207204 (2010) 
[24] S. Biermann et al., Phys. Rev. Lett. 94, 026404 (2005). 
[25] J. M. Tomczak and S. Biermann, Europhys. Lett. 86 3 (2009) 

37004, Phys. Rev. B 80, 085117 (2009). 
[26] Z. P. Yin, K. Haule and G. Kotliar, Nature Materials 10, 932 

(201 1), A. Kutepov et al., Phys. Rev. B 82, 045105 (2010). 
[27] K. Haule, J. H. Shim and G. Kotliar, Phys. Rev. Lett. 100, 

226402 (2008). 

[28] V. I. Anisimov et al., J. Phys. : Cond. Matter 21, 075602 (2009). 
[29] L. Craco et al., Phys. Rev. B 78, 134511 (2008). 
[30] D. Wu et al., Phys. Rev. B 73, 235206 (2006). 
[31] P. Werner, et al, Nature Physics 8, 331337 (2012). 



